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engines.  The  models  are  then  used  to  study  the  effect  of  changes 
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developed  by  Koenigsberg/[-2T^  - 

The  steady-state  behavior  of  both  models  is  essentially  the  same, 
with  the  analytical  model  consistently  more  conservative  than  the 
simulation  model. 

The  analytical  model  has  been  adapted  to  a  time-shared 
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immediately  obtain  effects  on  the  system  of  changes  in  any  one 
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I.  INTRODUCTION 

Litton  Industries  has  been  awarded  a  contract  to  build 
30  SPRUANCE  Class  destroyers  for  the  United  States  Navy  with 
first  delivery  scheduled  for  December  197^.  Each  ship  will 
be  equipped  with  four  General  Electric  LM-2500  gas  turbine 
engines  for  main  propulsion.  The  present  maintenance  plan 
provides  for  an  age  replacement  policy,  i.e.,  each  engine 
will  be  removed  upon  attaining  some  prescribed  number  of 
operating  hours  or  upon  failure,  whichever  occurs  first. 

Upon  removal,  the  engine  will  be  replaced  by  a  ready-for- 
issue  (RFI)  engine  obtained  from  an  RFI  pool  and  the  removed 
engine  sent  to  a  central  repair  facility  for  repair  and 
return  to  the  RFI  state,  thus  forming  a  closed-loop  use/- 
repair  system.  The  number  of  spare  engines  to  buy,  and  the 
repair  capability  to  establish,  are  to  be  determined. 

This  thesis  describes  and  compares  two  models  which  may 
be  used  as  tools  in  making  logistic  support  decisions  for 
closed-loop  systems,  such  as  the  use/repair  cycle  of  the 
LM-2500  gas  turbine  engines.  The  models  are  then  used  to 
study  the  effect  of  changes  in  various  system  parameters. 

The  models  are  a  computer  simulation  developed  by  Cushing, 
Gautier,  and  Long  [1]  and  an  analytical  model  developed  by 
Koenigsberg  [2]. 

The  main  decision  variables  are  the  pool  size,  number  of 
repair  lines,  time  to  repair,  and  the  ratio  of  time-between- 
overhaul  (TBO)  to  annual  hours  per  engine  (TOR). 

e 


With  TBO  =  4,000  hours,  TOR  *  2,000  annual  hours  per 
engine,  pool  size  =  19,  repair  time  =  5  weeks,  number  of 
repair  lines  =  6,  as  the  base  case  for  analysis,  the  para¬ 
meters  were  varied  through  ranges  abo^e  and  below  their  base 
case  values.  The  results  consistently  reveal  that  the  pool 
size  can  vary  over  a  wide  range,  from  a  minimum  of  12,  with 
little  or  no  effect  on  the  status  of  the  operating  fleet. 
However,  the  status  of  the  fleet  is  extremely  sensitive  to 
the  other  parameters,  especially  the  number  of  repair  lines, 
repair  time,  and  TBO.  For  example,  if  all  parameters  except 
TBO  are  kept  at  base  case  values,  and  TBO  is  increased  from 
4,000  to  6,000  hours,  the  average  number  of  engines  short 
in  the  fleet  drops  from  about  10  to  virtually  zero. 

Both  models  give  essentially  the  same  steady-state 
results  with  the  analytical  model  requiring  much  less  com¬ 
puter  time.  This  makes  possible  the  use  of  the  analytical 
model  in  a  real  time  interactive  computer  system.  However, 
the  simulation  model  provides  the  flexibility  to  handle  a 
larger  range  of  problems  such  as  those  with  transportation 
delays.  Because  of  the  long  running  time  and  other  inherent 
characteristics  of  the  simulation  model,  it  cannot  be  used 
in  a  real  time  mode. 

The  analytical  model  has  been  adapted  to  a  time-shared 
interactive  computer  system  which  allows  the  policy  planner 
to  immediately  obtain  effects  on  the  system  of  changes  in 
any  one  of  seven  parameters. 


This  thesis  contains  six  sections,  of  which  this  is 
the  first,  and  2  Appendices.  Section  II  contains  a  des¬ 
cription  of  the  system  and  a  description  of  the  simulation 
and  analytic  models.  Section  III  contains  a  comparison  of 
the  two  models  under  steady-state  conditions.  Section  IV 
contains  the  analysis  of  the  steady-state  results  obtained 
with  the  analytic  model  and  effects  of  transportation  de¬ 
lays  and  spare  engine  buying  plans  obtained  with  the  simu¬ 
lation  model.  Section  V  contains  the  explanation  of  the 
interactive  computer  model.  Section  VI  summarizes  this 
study.  Appendix  A  contains  sample  inputs  and  outputs  of 
the  interactive  computer  program.  Appendix  B  contains  the 
simulation  and  analytic  model  computer  programs. 


II.  THE  MODELS 


A.  SYSTEM  TO  BE  MODELED 

Before  starting  any  discussion  of  the  two  m<JRels,  we 
describe  the  system  to  be  modeled.  T^e"* SPRUANCE  CLASS 
destroyers  are  to  be  powered  by  the  Ge*fcj»l  Electric  LM- 
2501  gas  turbine  engines.  The  majnte'j^pce  policy  for  these 
engines  states  that  only  a  minimum  amount  of  onboard  main¬ 
tenance  is  to  be  performed  by  ship's  personnel.  Instead, 
the  engines  ^rkr'designed  for  easy  removal  and  replacement, 
so  that/  failed  engines  can  be  shipped  to  a  central  rework 
facility  for  maintenance.  This  sort  of  maintenance  policy 
requires  that  spare  engines  be  purchased  to  replace  those 
that  are  undergoing  rework.  The  purpose  of  this  is  to  keep 
the  fleet  operational  as:  much  of  the  time  as  possible  by 
replacing  failed  engines. 

The  typical  cycle  of  an  engine  begins  with  an  engine 
operating  in  the  fleet.  Based  on  a  policy  of  age  replace¬ 
ment,  engines  are  removed  from  the  fleet  either  when  they 
have  accumulated  a  specified  number  of  operating  hours,  or 
when  an  actual  failure  occurs  prior  to  the  specified  remov¬ 
al  time.  When  an  engine  has  been  removed,  it  is  transported 
to  the  rework  facility.  This  shipping  causes  a  delay  in 
the  engine's  arrival  at  the  rework  facility.  Once  at  the 
rework  facility  an  engine  may  or  may  not  begin  rework 
depending  on  whether  the  rework  facility  is  congested. 
However,  once  inducted  into  a  rework  line,  each  engine 


requires  a  known  amount  of  time  to  be  reworked.  When  rework 
has  been  completed,  the  engine  is  placed  in  a  ready-for- 
issue  status  to  await  a  requirement  for  use  in  the  fleet. 
There  is  usually  no  delay  in  transportation  from  the  reaay- 
for-issue  pool  to  a  ship,  since  a  ship  can  predict  in  advance 
when  an  engine  will  reach  its  scheduled  removal  time.  How¬ 
ever,  if  a  failure  occurs  early,  the  engine  would  experience 
a  delay  being  sent  to  the  ship.  A  block  diagram  of  this 
system  is  shown  in  Figure  1. 

B.  THE  SIMULATION  MODEL 

The  computer  simulation  model  is  a  General  Purpose 
Simulation  System/360  event  stored  routine.  This  model 
begins  operation  by  generating  a  transaction,  in  this  case 
the  first  ship  delivery,  which  sets  everything  into  motion. 
The  effects  of  this  transaction  are  as  follows: 

1.  It  causes  the  hotatable  pool  to  be  loaded  in  accor¬ 
dance  with  the  loading  instructions  provided  in  the  model. 
This  model  provides  for  loading  the  entire  quantity  at  time 
zero  or  at  any  other  time  during  the  simulation,  or  to  load 
increments  up  to  some  total  quantity  at  various  times  during 
the  simulation. 

2.  It  causes  the  simulation  clock  to  start.  The  model 
is  run  with  increments  of  one  week  as  the  basic  unit  of 
time,  i.e.,  the  parameters  are  all  scaled  to  weeks.  The 
system  is  observed  at  the  end  of  each  week  and  the  status  of 
the  operating  fleet,  rotatable  pool,  repair  lines,  etc.  is 
recorded. 
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3.  It  causes  operating  units  to  be  introduced.  The 
model  has  the  flexibility  to  provide  for  a  build-up  in  the 
number  of  units  operating  by  some  prescribed  delivery 
schedule,  e.g.,  introduction  of  four  engines  with  each 
DD-963  class  destroyer  in  accordance  with  the  ship’s  deliv¬ 
ery  schedule . 

Various  maintenance  plans  may  be  used,  from  the  pres¬ 
ently  planned  age-replacement  policy,  to  one  of  removal 
based  on  some  form  of  condition  monitoring.  In  each  case, 
a  function  determining  the  time  to  removal  of  each  engine 
must  be  provided  to  the  model. 

A  time  to  removal  is  assigned  to  each  engine  as  it 
arrives  in  the  operating  fleet  by  a  Monte  Carlo  technique, 
whether  from  new  ship  delivery  or  as  a  replacement  engine 
drawn  from  the  rotatable  pool. 

As  the  time  period  arrives  for  an  engine  to  be  removed, 
it  leaves  the  fleet  and  a  replacement  is  drawn  from  the 
pool,  if  available.  Otherwise,  replacement  will  be  made  with 
the  first  engine  available  at  the  pool.  The  removed  engine 
proceeds  to  the  repair  facility,  either  directly,  or  through 
a  transportation  delay,  if  required. 

Upon  arriving  at  the  repair  facility,  the  engine  is 
inducted  on  a  first-come-first-served  basis.  If  the  repair 
lines  are  all  busy,  the  engine  is  held  in  a  queue  until  a 
repair  line  becomes  available. 

The  engines  proceed  through  repair  in  accordance  with 
the  repair  function  which  assigns  a  time  to  repair  for  each 
engine  inducted. 
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Upon  completion  of  repair,  the  ready-for-issue  engines 
then  proceed  directly  to  the  rotatable  pool  where  they 
remain  until  called  for  by  the  requirements  of  the  operating 
fleet . 

The  parameters  which  must  be  specified  as  inputs  to  the 
simulation  are: 

1.  The  ship  delivery  schedule. 

2.  The  number  of  engines  per  ship. 

3.  The  function  to  generate  engine  removal  times. 

4.  The  rotatable  pool  loading  plan,  in  time  increments 
and  quantity  to  be  loaded  each  increment  of  time. 

5.  The  transportation  delays,  if  required. 

•f.  The  number  of  repair  lines. 

7.  The  function  describing  the  time  to  repair. 

8.  The  time  for  which  the  simulation  is  to  run. 

9.  The  standard  GPSS  statistics  desired  as  output,  or 
any  special  outputs  desired. 

The  model  typically  requires  about  2.5  minutes  of  com¬ 
puter  time  to  provide  statistics  for  ten  iterations,  each 
of  20  years  duration  for  parameters  in  the  ranges  of 
interest. 

C.  ANALYTICAL  MODEL 

The  LM-2500  gas  turbine  engine  use/rework  cycle  can 
also  be  described  in  terms  of  a  finite  cyclic  queue  model. 

We  make  the  following  three  assumptions: 
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1.  There  is  no  delay  in  the  transportation  of  engines  to 
and  from  the  rework  facility. 

The  repair  time  at  the  rework  facility  is  exponentially 
distributed  with  rate 
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where  T^  is  the  average  time  to  rework  an  engine. 

3.  The  time  an  engine  operates  in  the  fleet  is  exponen¬ 
tially  distributed  with  rate 
u2  *  1/T2, 

where  Tg  is  the  mean  number  of  weeks  required  for  an 
engine  to  accumulate  a  specified  number  of  operating 
hours . 

The  consequences  of  all  three  of  the  above  assumptions  will 
be  discussed  when  this  model  and  the  simulation  model  are 
compared  in  section  III.  Under  these  assumptions,  one  can 
determine  equations  for  the  various  system  performance 
measures . 

In  order  to  discurs  the  model  more  easily,  we  define  the 
following  parameters : 

N  *  The  total  number  of  engines  in  the  inventory.  This 
includes  both  the  number  of  engines  to  be  installed 
in  the  fleet  and  the  number  of  spare  engines. 

A  =  The  number  of  engines  to  be  installed  in  the  fleet. 

M  =  The  number  of  repair  lines  in  the  rework  facility. 
Let  us  now  define  two  stages  for  the  system.  An  engine 
must  be  in  one  of  these  two  stages.  Stage  1  (The  rework 
facility)  is  an  M-server  queue  each  with  an  exponential 
service  time  with  mean  rate  y-^.  When  all  M  repair  lines 
are  busy  an  incoming  engine  waits  in  a  queue  for  rework. 
Stage  2  (The  fleet)  is  an  A-server  queue  each  with  an  ex¬ 
ponential  service  rate  y,,.  When  all  ships  have  their  full 
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complement  of  engines,  a  ready-for-issue  (RFI)  engine  waits 
in  a  pool  until  there  is  a  demand  from  the  fleet  for  an 
engine . 

.  The  two  stages  form  a  closed  system  in  that  no  engines 
enter  or  leave  the  system.  Since  N  will  remain  constant, 
we  can  define 

P.^  *  Probability  there  are  i  engines  in  stage  1  (hence, 
N-i  engines  in  stage  2). 

For  the  finite  cyclic  queue  just  described,  where  the 
number  of  spare  engines  is  greater  than  or  equal  to  the 
number  of  repair  lines  in  the  rework  facility,  i.e. 

N-A  <  M, 

Koenigsberg  in  reference  [2]  determines  the  P.^  to  be: 

P±  =  (Ai/i!)(p2/v1)1  PQ,  0  <  i  <  M,  (la) 

=  (A1/!*!  M1^)  (pg/p^1  PQ,  M<i<N-A,  (lb) 

«  (A!AN“A/(N-i)!M!Mi"M)(p2/p1)i,  N-A<i<N.(lc) 

It  is  also  necessary  that 
N 

E  P,  «  1.  (2) 

i»0  x 

These  P^  will  be  used  to  calculate  the  various  system 
measures. 

Writing  the  ith  equation  as 


<« -A- ■  ^ •  »E -4^9939^  S7f*, - 


^^vFj£F‘v»  -», «  - 


«. 4^'a»  ■$'*•*  tri***  »« 


then  summing  over  i  on  both  sides  of  (3)  gives 
N 

P0  z  at  -  1.  (4) 

0  i-0  1 


Therefore 

N 

?0  -  1/  £  a,.  (5) 

0  i-0  1 


Substituting  (5)  into  (3)  we  obtain 

N 

P,  »  a,/  £  a.,  i  «  0,1,..., N.  (6) 

1  1  i=0  1 


Direct  computation  of  the  a^’s  leads  to  numerical 
problems.  For  example,  let  N  *  1^0,  M  *  6,  and  A  »  120. 
Then 

_  1201  (120)20 
120  20!  61  e11"  ' 

Direct  computation  of  such  a  number  in  any  machine  will 
lead  to  problems  since  100!  is  of  the  order  of  lO1-^.  Thus 
we  use  the  following  recursive  forms  to  calculate  these 
coefficients : 

Note  that  aQ  =  1,  and  let  p  *  (yg/y^. 

Then 

a.  ,  A  p 

ai  =  - -  ,  1  <  i  <  M,  (7a) 

a.  ,  A  p 

-  -±=i -  ,  M  <  i  <  N-A,  (7b) 
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a1_1  p  (N+2-i) 


M 


N-A  <  1  <  N. 


(7e) 


In  reference  [2],  Koenigsberg  defines  certain  measures 
of  system  performance.  The  measures  are  interpreted  in 
terms  of  the  LM-2500,  use/rework  cycle. 

p^  «  Mean  number  of  engines  in  the  rework  facility  ■ 
and  adjoining  waiting  line. 


N 


z  i  P. . 
i-G  X 


ng  *  Mean  number  of  engines  in  the  fleet  and  adjoining 
RFI  pool. 


-  N  -  nr 


w,  »  Mean  number  of  engines  awaiting  rework. 


N 

z 


i=M+l 


(i-M)  P±. 


w2  =  Mean  number  of  RFI  engines. 


N-A 


*  Z  (N-A-i)  P,  . 
i-0  x 


d^  -  Percent  time  all  repair  lines  are  busy, 


N 


-  (E  P. )  x  100 
i-M  x 


*  Mean  number  of  engines  being  repaired. 


"  nl  -  W1  * 
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H2  ■  Mean  number  of  engines  operating  in  the  fleet, 

-  n2  -  w2  . 

Um  *  Mean  percent  rework  capacity  used, 

»  (Hx/M)  x  100. 

Un  *  Mean  percent  of  total  engines  in  operation  in 
the  fleet, 

=  (H2/N)  x  100. 

Three  additional  system  measures  not  given  by  Koenigs- 
berg  in  reference  [2]  are  also  included  in  the  output.  The 
first  is  cycle  time  (C).  This  is  the  average  time  required 
for  an  engine  to  complete  a  full  cycle  of  the  system.  To 
derive  this  measure,  we  first  note  that  the  rate  into  a 
stage  is  equal  to  the  rate  out  of  the  same  stage.  Then  in 
steady  state  the  rate  out  of  the  fleet  stage  is  y2  Hg.  This 
is  the  rate  everywhere  in  the  system,  so  the  cycle  time  (C) 
is  given  by 

C  =  n/v2h2  . 

To  calculate  the  probability  the  fleet  is  short  i  or  more 

engines  (P_(i)),  we  use  the  following: 
s 

N-A+i 

P_(i)  =  1  -  l  P,. 

8  i=0  1 
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To  calculate  the  average  number  of  engines  short  in  the 
fleet  (?) ,  we  use 


s  *  A  -  H 


2  * 


All  of  the  above  measures  of  system  performance  are  printed 
in  the  output  exeept  n^  and  n2  which  are  only  used  to  cal¬ 
culate  other  measures. 

v  There  are  seven  parameters  of  the  LM-2500  use/rework 
cycle  which  serve  as  the  inputs  to  the  model.  They  are: 

1.  Number  of  ships  (S), 


2.  Number  of  engines  per  ship  (E), 

3.  Number  of  spare  engines  (P), 

*1.  Number  of  repair  lines  (M) , 

5.  Average  rework  time  per  engine  in  weeks  (T^), 

6.  Time  between  overhauls  (TBO), 

7.  Engine  operating  hours  per  year  (TOR). 


The  above  seven  parameters  are  converted  internally  in 
the  program  to  the  actual  parameters  used  by  the  model  in 
the  following  way: 


A  c  S  x  E, 
N  *  A  +  P, 
V1  *  1/TX, 


T2  =  (TBO/TOR)  x  52  x  0.9, 


1/T, 


M  is  unchanged. 

All  of  the  above  are  self  explanatory  except  T? .  This 
equation  for  T2  is  explained  fully  in  the  section  on  compar 
ison  of  the  models.  Briefly,  it  calculates  a  mean  time  to 
replacement  consistant  with  the  age  replacement  policy  for 
engine  removals. 
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III.  COMPARISON  OF  MODELS 


The  simulation  model  was  developed  because  ”...  analysis 
of  the  DD-963  inventory  problem  disclosed  that  system  com¬ 
plexities  rendered  analytical  modeling  infeasible.”  (See 
[1]  ,  page  17)  That  is  true  if  all  of  the  transient  prob^ 
lems  are  to  be  included.  However,  inclusion  of  the  rela¬ 
tively  short  transient  period  in  a  long-term  program  could 
lead  to  decisions  which  are  harmful  in  the  long-term.  For 
that  reason,  steady-state  results  from  the  simulation  and 
the  analytical  model  were  compared  to  determine  any  differ¬ 
ences  in  results  from  the  two  models.  Since  the  analytical 
model  has  certain  restrictions,  e.g.,  no  capability  to 
handle  transportation  delays,  the  simulation  model  was 
modified  to  eliminate  transportation  delays  and  to  elimin¬ 
ate  the  transient  effects  by  introducing  the  30  ships  and 
loading  the  pool  at  time  zero. 

The  two  models  were  run  with  the  same  parameters  over 
various  ranges.  The  results  of  the  two  models  were  extremely 
close,  with  the  analytical  model  consistently  more  conser¬ 
vative  than  the  simulation  model.  Most  important,  as  shown 
by  Figures  2,  3,  and  4,  the  response  to  changes  in  particu¬ 
lar  parameters  followed  the  same  pattern  in  the  two  models. 

As  an  example,  the  average  number  of  engines  awaiting  repair 
plotted  against  rotatable  pool  size,  with  the  4,000  hour  TBO, 
2,000  annual  hours  per  engine,  repair  time  of  5  weeks,  and 


6,  and  8  repair  lines,  is  displayed  in  Figure  2.  While 
the  number  waiting  is  slightly  less  with  the  simulation 
model  than  with  the  analytical  model,  both  models  show  the 
same  behavior  with  a  change  in  pool  size.  For  a  repair  fac¬ 
ility  with  only  4  repair  lines  and  a  pool  size  of  15,  the 
waiting  line  is  long,  42  for  the  simulation  model  and  46 
for  the  analytical  model.  However,  more  important  is  that 
each  additional  engine  added  to  the  pool  increases  the  size 
of  the  waiting  line  by  one  engine  for  both  models.  Similar 
results  are  evident  with  the  other  comparisons  in  Figures 
3  and  4. 

Although  from  the  foregoing  model  comparison  we  conclude 
that  both  give  the  same  steady-state  results,  it  should  be 
pointed  out  that  the  time  in  the  fleet  and  repair  time 
distributions  used  in  the  two  models  were  very  different. 
This  shows  that  the  analytical  mode  is  quite  robust  and 
deviation  from  the  exponential  distributions  do  not  lead 
to  serious  errors. 

In  the  simulation  model,  the  time  an  engine  stays  in 
the  fleet  has  an  empirically  determined  distribution  func¬ 
tion,  constructed  by  Cushing,  Gautier,  and  Long  [1]  using 
actual  records  of  DD  type  operations.  The  analytical  mode] 
has  an  exponential  distribution  function  with  a  mean  time 
in  the  fleet  determined  by  the  TBO  and  TOR.  The  repair 
time  for  the  simulation  is  deterministic,  whereas  the 
analytical  model  uses  an  exponential  distribution  function 
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with  a  mean  repair  time  equal  to  the  deterministic  value 
ol*  .the  simulation.  1  > 

In  our  comparison  of  the  distributions  of  the  time  an 
engine  stays  in  the  fleet,1  the  4, ,000  hour  TBO  and  2,000  TOR 

>'  i 

operating  profile  are  used.  These  result  in  a  mean  time 

*  •  i  i 

i  in  the  fleet  of  106  weeks  which  is  used  in  the  exponential 
distribution. 

i 

As  previously  explained  in  the  simulation  model,  a 

»  *  • 

!  *  S  I  ,  * 

time,  in  weeks,  for  an  engine  to  react!  4,000  hours  is  a 


I 


! 


■random  variable  which  we, will  define  to  be  t.  The  distri- 

1 

bution  of  X  is  discrete  with  the  following  mass  function: 


1  x 


66  weeks 
72 

78 

ti 

84 

1  90 

’ll 

96 

11 

102 

"  . 

108' 

n 

114 

it 

1|20 

i« 

126 

ti 

132 

it 

138, 

11  '■ 

144 

11 

150 

1  ti 

156 

it 

162 

ii  ' 

168 

■  it 

174 

it 

elsewhere 


'  i 

Once  X  is  determined  for  an 

i 

probability  that  the  engine  will 


P(X=X) 


.009 

.019 

.019 

.027  1 

.046  1 
.055 
.076 
.093 
.120 
.120 
>.120 
.102 
.074 
.037 
.029 
.019 
.019  , 

.009 

.009 

.000 

engine,  ithere  is  a  0.8 

I 

,  be  removed  at  that  time  and 


a i  0.2  probability  that  it  'will  .fail  before  the  time  X.  .The 
1 

1  '  *  1  ’  •  1 
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distribution  of  the  time  for  the  early  failure  is  assumed 
to  be  uniform  between  zero  and  X  weeks. 


The  distribution  of  the  time  in  the  fleet  is  a  random 
variable  Y,  a  function  of  the  random  variable  X,  the  U(0,X) 
random  variable,  and  the  probability  of  reaching  the  de¬ 
sired  TBO.  The  conditional  distribution  of  Y  given  X  is: 


P(Y<y|x=x)  » 


for  y  <  x, 
for  y  >  x. 


Figure  5  is  a  sample  of  the  conditional  distribution  of  Y. 


Figure  5.  Sample  Conditional  Distribution  of  Y  given  X 


The  unconditional  distribution  of  Y  is 

y  174 

P(Y  <  y)  =  E  p(x)  +  E  p(x)  . 

x=0  x=y+l  x 

A  plot  of  this  distribution  function  is  displayed  in 
Figure  6  along  with  the  exponential  distribution  function 
with  the  same  mean. 

While  the  two  functions  are  clearly  quite  different, 
the  probability  of  a  removal  prior  to  128  weeks  of  operation 
is  approximately  0.7*1  for  both  distributions.  There  is  a 
probability  of  approximately  .16  that  an  engine  stays  in 
the  fleet  in  excess  of  17*1  weeks  using  the  exponential 
distribution  which  is  zero  for  the  empirical  distribution. 
This  provision  for  times  greater  than  17*1  weeks  prior  to 
removal  for  some  engines  has  a  "real  world"  appeal. 

This  comparison  was  done  with  the  *1,000  hour  TBO  and 
the  2,000  annual  hours  per  engine  TOR.  The  mean  time  in 
the  fleet  derived  from  the  empirical  distribution  function 
for  this  combination  is  106  weeks.  This  is  derived  by 
finding  E[X]  and  then  finding  E[Y]  as  follows: 

E[Y]  =  E[E[Y|X]]  =  EC'^P*  +  .8X3  =  .9  E[X]  . 

Cushing,  Gautier,  and  Long  [1]  developed  empirical 
distribution  functions  for  various  TBO  and  TOR  combinations. 
Rather  than  calculate  the  mean  of  each  of  these  functions 
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for  use  in  the  analytical  model,  we  have  defined  the 
following : 

E[X]  (in  weeks)  “  52(TBO/TOR)  and 

E[Y]  (in  weeks)  =  .9E[X]  ,  for  all  TBO  and 

TOR  combinations. 


If  the  interest  is  in  the  number  of  engines  that  have 
to  be  removed  at  exactly  30  weeks  of  operation  in  the 
fleet,  or  some  other  such  specific  time,  the  two  models 
would  not  necessarily  agree.  But  if  the  primary  interest 
is  in  analyzing  overall  results,  e.g.,  the  average  number 
of  engines  short  in  the  fleet,  then  the  two  models  give 
virtually  the  same  results. 

Given  the  current  state  of  the  system  ,  no  one  knows 
which  distribution  will  best  describe  the  actual  operations. 
But  for  analyzing  the  long  term  behavior  of  the  whole 
system,  either  model  could  be  used.  The  choice  of  which 
model  to  use  should  be  determined  by  the  advantages  one 
might  have  over  the  other. 

The  analytical  model  is  limited  in  that  it  provides  only 
steady-state  results.  For  this  reason,  the  simulation  is 
used  to  analyze  the  transient  time  problems  of  incremental 
ship  delivery,  incremental  pool  purchases,  and  the  addi¬ 
tional  problem  of  transportation  delays.  The  analytical 
model  provides  a  rapid  efficient  method  of  measuring  the 
effects  of  changing  various  policy  variables.  The  results 
obtained  on  the  computer,  using  the  analytical  model. 
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require  about  2  seconds  computer  time  for  each  set  of 
parameters  investigated,  whereas  the  simulation  requires 
about  2.5  minutes  tor  the  same  analysis. 
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IV.  INTERACTIVE  COMPUTER  PROGRAM 

To  provide  the  policy  planner  with  a  real-time  system 
for  determining  the  effects  of  various  parameter  changes 
on  the  LM-2500  use/repair  cycle,  an  interactive  computer 
program  has  been  developed  which  uses  the  analytical  model. 

To  operate  the  interactive  computer  system,  the  policy 
planner  would  require  the  use  of  a  remote  terminal  which 
can  be  connected  to  the  computer  in  which  the  program  is 
stored.  He  would  then  go  through  a  standard  procedure  i;o 
begin  the  execution  of  the  program.  This  procedure  is  a 
function  of  the  computer  center  so  it  will  not  be  explained 
here. 

Once  execution  begins,  the  program  will  print  instruc¬ 
tions  for  entering  the  use/repair  cycle  parameters.  Seven 
parameters  are  to  be  entered.  They  are: 

1.  Number  of  ships 

2.  Number  of  engines  per  ship 

3.  Number  of  spare  engines 

4.  Number  of  repair  lines  in  the  rework  facility 

5.  Average  rework  time  per  engine  in  weeks 

6.  Time  between  overhauls  in  hours 

7.  Engine  operating  hours  per  year. 

When'  all  of  these  parameters  have  been  entered,  the  program 
will  type  out  the  parameters  and  their  value  so  that  these 
can  be  checked  for  possible. errors.  If  they  are  in  fact 
correct,  the  user  types  a  1  (one),  and  the  program  will 


$ 

I 


£ 

f 


calculate  the  use/repalr  cycle  system  performance  measures. 
If  they  are  not  correct  the  user  can  type  a  0  (zero),  and 
the  program  will  return  to  lt^-s  start  so  the  correct  param¬ 
eters  can  be  entered.  In  Appendix  A,  on  page  50,  an  actual 
example  of  this  Input  procedure  is  included. 

On  page  52,  a  sample  of  the  output  is  also  listed. 

This  whole  process  of  Input  and  output  requires  about  five 
minutes.  The  decision  maker  can,  if  he  chooses,  restart 
the  program  with  a  new  set  of  parameters.  Actual  computer 
time  per  set  of  parameters  is  less  than  two  seconds. 
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V.  ANALYSIS 

A.  STEADY  STATE  RESULTS 

Of  the  many  measures  of  system  performance,  the  follow¬ 
ing  measures  were  chosen  as  providing  the  most  information 
to  the  policy  planner: 

1.  Average  number  of  engines  short  in  the  fleet  (s’), 

2.  Average  number  of  RPI  engines  (w2), 

3.  Average  percent  of  rework  capacity  used  (Um), 

4.  Average  number  of  engines  awaiting  rework  (w^). 

When  looking  at  these  measures  for  various  combinations  of 
system  parameters,  it  must  be  kept  in  mind  that  they  are 
steady  state  results.  That  is,  this  is  how  the  system  will 
look  after  It  has  been  operating  over  a  long  period  of  time. 
The  effects  of  the  start-up  of  the  system  have  little  or 
no  effect  on  the  value  of  the  steady  state  system  perform¬ 
ance  measures. 

In  figures  7  through  10,  the  four  measures  above  are 
plotted  versus  the  ratio  of  TBO  to  TOR.  All  four  graphs 
are  based  on  19  spare  engines.  The  effects  of  a  change  in 
the  number  of  spare  engines  are  investigated  below.  Table  1 
converts  a  given  TBO  and  TOR  into  the  ratio  shown  on  the 
graphs.  In  each  figure,  a  set  of  nine  curves  is  plotted. 
Each  curve  corresponds  to  a  given  number  of  repair  lines 
(M)  and  an  average  repair  time  (T^.  Each  line  is  labeled 
(M,TX). 
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TABLE  1 


Converts  a  given  TBO  and  TOR  to  tie  ratio  used  in  figures  7 

through  10 


TBO  TOR 


1650 

2000 

2500 

4000 

2.4 

2.0 

1.6 

5000 

3.0 

2.. 

2.0 

6000 

3.6 

3.0 

2.4 

7000 

4.2 

3.5 

2.8 

8000 

4.8 

4.0 

3.2 

9000 

5.5 

4.5 

3.6 

10000 

6.1 

5.0 

4.0 

11000 

6.7 

5.5 

4.4 

12000 

7.3 

6.0 

4.8 
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Figure  7  shows  the  relation  between  average  number  of 
engines  short  in  the  fleet : (?)  and  ^BO/TOR.  For  each  of 
the  nine  curves  there  is  a  dramatic  decrease  in  s’  as  TBO/TOR 
.increases  to  about  4,  usually  an  inbrease  in  TBO/TOR  is 

|  i 

brought  aDout*by  An  increase  in  engine  reliability.  Also 

»  *  i 

, it  can  be  seen  that  there  is  a  strong  effect,  due  to  the 
combination  of  M  and  T^.  That  is,  when  there  is  adequate 

■  i 

rbpair  capacity  in  terms  of' size  and  speed,  then  the  num¬ 
ber  of  engines  short  in  the  fleet  is  reduced  rapidly. 

:  I  . 

The  plots . of  the  other  three  measures  in  figures  8,  9* 
and  10  show  the  same  dramatic  changes  as  TBO/TOR  is  in- 

!  i 

.  ,  1  •  . 

creased  to  about  4.  The  combination  of  M  and  T,  also  pro<- 

•  !  1  '  X 

.  duce  the  same  sharp  effects  as  the  rework  facility  is  made 

more  Oapable  of  reworking  engines. 

!  •  1 

■  For  example,,  let  us  look  at  a  system  where 

i  M  *=  6, 

!  :  ,  •  •  1 

A1  i 

•  TBO  =!  H000, 

i  !  ,  i 

TOR  *  2000. 

•  ! 

We  will  then  change  TBO  to  6000  and  look  at  the  result  -, ,  of 
this  change  in  TBO/TOR  from  2  to  3*  ,Now  enter  each  figure 
with  these  two  ratios.  'In  figure:  7,  entering  at  TBO/TOR 

*  t  t 

=  2  >we  find  the  value  of  ?  to  be  10  on  the  curve  labeled 

i  »  i  » 

I  S 

.(6,5)  •  When  entering  the  graph  at  TBO/TOR  =  3>  we:  find 
that  'the  (6,5)  curve  has  gone  to  0  for  a,  ratio  less  than- 3. 
So  the  change  in  TBO  pas  resulted  in  a  reduction  of  the 
'  number  Of  engines  short  in  the  fleet  by  10.  Entering 
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figure  8  in -the  same  manner,  shows  that  the  number  of  RFI 
engines  is  increased  from  2  to  about  13.  Prom  figure  9> 
we  see  that  the  average  facility  capacity  used  decreased 
from  99 %  to  73Jt  -  Prom  figure  10,  we  see  that  the  average 
number  of  engines  awaiting  rework  decreased  from  20  to  less 
than  1.  It  is  therefore  quite  evident  that  any  program  to 
increase  the  reliability  of  the  engines  will  be  of  great 
benefit  to  the  operating  fleet  and  will  cut  down  on  the 
required  rework  capacity. 

The  above  graphs  were  produced  using  19  spare  engines, 
with  a  fleet  of  120  engines.  A  natural  question  is  what 
effect  will  a  change  in  the  number  of  spare  engines  have 
on  the  four  system  measures?  For  this  reason,  the  four 
measures  are  plotted  versus  the  number  of  spare  engines  in 
figures  11  through  14.  Each  figure  uses  a  value  of  T2  -  106 
weeks,  M  =  6,  and  T^  =  5  weeks.  In  some  figures,  other 
values  of  M  and  T.^  are  used  and  will  be  explained  as  they 
are  needed. 

Figure  11  shows  the  average  percent  of  facility  capacity 
used  (Um)  versus  spare  engines.  This  measure  is  insensitive 
to  changes  in  the  number  of  spare  engines.  So  when  com¬ 
paring  Um,  the  number  of  spare  engines  need  not  be  consid¬ 
ered. 

In  figure  12,  the  average  number  of  engines  short  in 
the  fleet  is  plotted.  This  measure  does  show  some  effects 
due  to  the  number  of  spare  engines.  However,  the  effect 
is  larger  between  8  and  15  spare  engines  than  between  15 
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and  23  spare  engines.  When  changing  the  number  of  spares 
from  8  to  15,  the  number  short  is  reduced  by  only  2.  In 
the  range  from  15  to  23,  this  reduction  is  only  one  engine. 
Therefore,  an  increase  in  the  number  of  spare  engines  has 
very  little  effect  in  reducing  shortages  in  the  fleet. 
Figure  7  shows  that  there  is  a  better  return  received  by 
developing  the  right  rework  facility  than  buying  spare 
engines . 

Figure  13  shows  a  different  sort  of  result  for  the 
average  number  of  RFI  engines.  Whereas  the  average  number 
of  engines  short  did  not  change  much,  this  measure  does. 

In  this  case,  three  curves  are  plotted,  the  standard  (6,5) 
case,  the  best  case  (8,4),  and  the  worst  case  (4,6).  The 
number  of  RFI  engines  is  low  and  unchanging  for  the  (4,6) 
case.  When  the  rework  facility  capacity  is  increased 
(increase  M,  decrease  T^),  the  number  of  RFI  engines  begins 
to  increase  with  the  number  of  spare  engines.  However  it 
is  the  parameters  of  the  rework  facility  that  are  causing 
the  effects,  not  the  number  of  spare  engines.  If  the  re¬ 
work  facility  is  large  enough  and  the  rework  time  short 
enough  then  the  spare  engines  end  up  in  the  RFI  pool. 

Figure  14  shows  that  the  waiting  line  at  the  facility 
does  increase  with  an  increase  in  the  number  of  spare 
engines.  This  is  perhaps  the  best  evidence  that  if  there 
is  not  enough  rework  capacity,  an  increase  in  the  number 
of  spare  engines  will  not  help  the  operating  fleet.  The 
engines  merely  add  to  the  length  of  the  rework  waiting 
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line.  The  slope  of  the  curves  approaches  1  as  the  rework 
facility  parameters  approach  their  worst  case.  So  in  the 
long  run,  to  buy  spares  to  keep  the  fleet  operating  will 
not  accomplish  its  purpose  unless  it  is  accompanied  by  a 
corresponding  increase  in  rework  facility  capacities. 


B.  TRANSPORTATION  DELAYS 

The  analytical  model  could  not  include  the  delays  which 
are  incurred  between  removal  from  the  ship  and  arrival  at 
the  repiir  facility,  and  the  similar  delay  incurred  between 
the  RFI  pool  and  the  fleet  when  an  unscheduled  removal  is 
experienced.  For  these  reasons,  the  simulation  model  was 
used  to  determine  the  effects  of  transportation  delays. 

The  simulation  model  was  run  with  the  120  engine  oper¬ 
ation,  a  pool  size  of  19,  repair  time  of  5  weeks,  6  repair 
lines,  and  transportation  delays  in  both  of  the  above  cases 
of  0,2,3,^,  and  8  weeks.  The  average  number  of  engines 
short  in  the  fleet  was  observed  in  each  case  with  the 


following  results : 


Ave.  #  engines  short  in  fleet 


0.3 

1.7 

2.8 
3.8 
8.3 


This  indicates  that  each  week  of  transportation  delay 
costs  the  operating  fleet  an  average  of  one  engine.  This 
might  be  interpreted  to  mean  that  an  increase  in  pool  size 
of  one  for  each  week  of  transportation  delay  would  counteract 


the  effect  of  1  week  transportation  delay.  To  test  this, 
for  the  transportation  delay  of  eight  weeks,  the  pool  size 
was  increased  from  19  to  27  and  the  simulation  repeated. 

The  average  number  of  engines  short  decreased  to  4.3  from 
8.3.  Thus,  two  additional  engines  in  the  pool  were  required 
to  reduce  the  number  short  in  the  fleet  by  1. 


C.  TRANSIENT  RESULTS 

Twc  different  spare-engine  buy  plans  were  tested  with 
the  simulation  model.  Each  plan  was  for  a  total  pool  size 
of  19.  The  plans  were  tested  with  the  4,000  hour  TBO, 

2,000  annual  operating  hours  per  engine  TOR,  repair  time  of 
5  weeks,  6  repair  lines,  transportation  delay  of  two  weeks, 
and  the  DD-963  planned  delivery  schedule.  The  buy  plans 
tested  were: 


Delivery  Dates 
Test  1  Test  2 


Number  of  engines 


Option  1  Option  2 


1  Jul  75  1  Jan  7 6 

..  76  »  77 

"  77  "  78 

.«  78  h  79 

"  79  "  80 


Using  the  percentage  of  time  there  are  no  engines  in 
the  pool  over  a  twenty  year  period  as  a  measure  to  compare 
the  effects  of  the  two  buying  plans,  the  following  results 
were  obtained: 


Option  1,  Test  1  92.75? 
Option  2,  Test  1  91.95? 
Option  1,  Test  2  92. 8j? 
Option  2,  Test  2  92. 0J? 


^^.rij^v-^r ;te> -**3-  ^^;*£r>ST^sa^wf ,x*?i^f«  4 .  «?v  ^£/*r  jyi^scpty; 


This  indicates  that  the  system  is  relatively  insensitive 


to  the  differences  in  these  two  buy  plans  over  reasonably 


long  time  periods.  The  particular  plan  selected  should 


probably  be  determined  by  economic-  considerations  in  the 


short  term,  say  the  first  5  years.  It  was  felt  that  the 


simulation  model  did  not  provide  reliable  analysis  of  the 


short  term  effects  of  the  two  buy  plans. 
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VI.  SUMMARY 

Using  as  measures  a)  the  average  number  of  engines 
short  in  the  fleet,  b)  the  average  number  of  engines  in  the 
RPI  pool,  c)  the  percentage  of  repair  capacity  used,  and 
d)  the  average  number  of  engines  awaiting  rework,  the 
significant  findings  of  this  study  can  be  summarized  as 
follows : 

1.  In  the  long  run,  the  four  measures  are  extremely 
sensitive  to  the  ratio  of  time  between  overhaul  to  annual 
engine  operating  hours.  The  system  is  improved  by  a  higher 
ratio. 

2.  In  the  long  run,  the  four  measures  are  extremely 
sensitive  to  the  repair  facility  capacity  as  measured  by 
the  number  of  repair  lines  and  average  engine  repair  time. 

3.  In  the  long  run,  the  percentage  of  repair  capacity 
used  and  the  average  number  of  engines  short  in  the  fleet 
are  relatively  insensitive  to  the  number  of  spare  engines 
purchased.  However,  the  average  number  of  engines  in  the 
RFI  pool  and  the  average  number  of  engines  awaiting  rework 
are  extremely  sensitive  to  the  number  of  spare  engines 
purchased.  The  additional  spares  end  up  in  an  Inactive 
status  over  the  long  run.  If  there  is  excess  repair  capac¬ 
ity,  they  end  up  in  the  RPI  pool.  If  there  is  insufficient 
repair  capacity,  they  end  up  in  the  repair  facility  waiting 
line.  For  cases  in  between  they  are  split  between  the  RFI 
pool  and  the  repair  facility  waiting  line. 


48 


^^5233 


i».  The  interactive  computer  system  developed  for  use 
with  the  analytical  model  provides  a  useful  method  for  the 
policy  planner  to  immediately  see  the  effects  of  changes 
in  system  parameters. 


-~>j.  „> 
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APPENDIX  A 

INTERACTIVE  COMPUTER  PROGRAM  INPUT/OUTPUT  FORMAT 
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APPENDIX  B 

COMPUTER  PROGRAMS 
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